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Abstract 

We consider a continuous time linear multi-inventory system with unknown de- 
mands bounded within ellipsoids and controls bounded within ellipsoids or poly- 
topes. We address the problem of e-stabilizing the inventory since this implies some 
reduction of the inventory costs. The main results are certain conditions under 
which e-stabilizability is possible through a saturated linear state feedback control. 
All the results are based on a Linear Matrix Inequalities (LMIs) approach and on 
some recent techniques for the modeling and analysis of polytopic systems with 
saturations. 
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1 Introduction 



We consider a continuous time linear multi-inventory system with unknown demands bounded 
within ellipsoids and controls bounded within ellipsoids or polytopes. The system is modelled 
as a first order one integrating the discrepancy between controls and demands at different 
sites (buffers). Thus, the state represents the buffer levels. We wish to study conditions under 
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which the state can be driven within an a-priori chosen target set through a saturated hnear 
state feedback control. Let £ be a maximal dimension of the target set, the above problem 
corresponds to £-stabilizing the state. 

Motivations for £-stabilizing the state derive from the benefits associated to keeping the state 
and consequently also the inventory costs bounded. This work is in line with some recent 
literature on robust optimization [1,5] and control [2] of inventory systems. Here as well as 
in [2] we focus on saturated linear state feedback controls since such controls arise naturally 
in any system with bounded controls. 

The main results of this work can be summarized as follows. Initially we introduce the 
necessary and sufficient conditions for the e-stabilizability in the form of an inclusion between 
convex sets. In the case where both demands and controls arc bounded within polytopes, it 
is well known that verifying such conditions is NP-hard [10]. Here, we prove that verification 
becomes easy when both demands and controls are bounded within ellipsoids (we will refer 
to it as the ellipsoidal case) . This is possible by rewriting the inclusion between ellipsoids in 
terms of unconstrained quadratic maximization. 

For the ellipsoidal case, we first characterize invariant sets through a fourth degree condition. 
As verifying such a condition is difficult, we then propose the best quadratic approximation 
of the same condition. We proceed by describing the region of linearity of the control and 
conclude by providing LMI conditions on the target set under which the saturated control 
^-stabilizes the system. The case where demands are bounded within ellipsoids and controls 
are bounded within polytopes (we will refer to it as the polytopic case) is an open problem 
and we propose certain sufficient LMI conditions to solve it. 

All the results are based on a Linear Matrix Inequalities (LMIs) approach in line with the 
recent work [6] on inventory/manufacturing systems. In particular, when addressing the poli- 
topic case, we use the same technique provided in [9] to rewrite the model with saturations 
in polytopic form. Once we do this, we can apply the LMI analysis covered in the book [7] 
for polytopic systems. 

This paper is arranged as follows. In Section 2, we formulate the problem. In Section 3, we 
introduce necessary and sufficient conditions for the admissibility of the problem. In Sections 4 
and 5 we study the problem with ellipsoidal and polytopic constraints respectively. Finally, 
in Section 6, we draw some conclusions. 



2 Problem Formulation 



Consider the continuous time linear multi-inventory system 

x{t) ^Bu{t) -w{t), 
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Fig. 1. Graph with one node and two arcs. 

where x{t) e IR" is a vector whose components are the buffer levels, u{t) e H"^ is the 
controlled flow vector, B e Q^x"^^ -^vith m > n and rank{B) — n is the controlled process 
matrix and w(t) G IR"^ is the unknown demand. To model backlog x{t) may be less than zero. 
Demands are bounded within ellipsoids, i.e., 

w{t) eW^{weW : w'^R^w < 1}. (2) 

In a first case, in the following referred as ellipsoidal case, controls are bounded within 
ellipsoids, 

u{t) eU ^{ueW^ : u^R^u < 1}. (3) 
In a second case, in the following referred as polytopic case, controls are bounded within 
polytopes 

u{t) e U = {u eR"" : u- < u < u+} (4) 

with assigned u~. 

For any positive definite matrix P e M"^", define the function V{x) — x'^Px and the 
ellipsoidal target set 11 = {x G IR" : V{x) < 1}. In addition, for any matrix K G R"^", 
define as saturated linear state feedback control any policy 

, -Kx iiKxeU 
u = -sat{Kx} = { (5) 
u{x) G dU otherwise 

where hereafter dF indicates the frontier of a given set F. 

Problem 1 (e- stabilizing) Consider a system (1) in the ellipsoidal or polytopic case. Find 
conditions on the positive definite matrix P G M"^", under which there exists a saturated 
linear state feedback control u — —sat{Kx} such that it is possible to drive the state x{t) 
within the target set 11. 

Solving the above problem corresponds to £-stabilizing the state x within 11. 

Example 1 Throughout this paper we consider, as illustrative example, the graph with one 
node and two arcs depicted in Fig. 1. The incidence matrix is B — [1 1]. The continuous 
time dynamics is 

\i{t) 

-W = Ui{t) + U2{t) - w{t). 



x{t) = [l 1] 



U2{t) 



with demand bounded in the ellipsoid 



< 1 
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and with the following either ellipsoidal or polytopic constraints on the control u 

{Ul + U2 f < 1, (6) 

- 2 < lii < 3, -2 < lis < 1. (7) 
Finally, the target set is the sphere of unitary radius U — {x & M. : < 1}. 

3 Stability necessary and sufficient conditions 

System (1) is £-stabilizable if and only if for all w e W, there exists u e int{U} such that 
Bu — w (see, e.g., [3]). For the short of notation, the previous condition is usually expressed 
as 

BU D W. (8) 

Deciding whether (8) holds is NP-hard, when U and W are polytopes. Here, we prove that 
verifying (8) becomes easy when both U and W are eUipsoids. Observe that we can rewrite 
Bu — w as ub — B~^w — B~^Nun, where B — [B\N] being B a basis of B and N the 
remaining columns of B, correspondingly ub are the n components of u associated to the 
basis B and mjv are the m — n components of u associated to the columns in N . 

As we observe that (8) is equivalent to 

max min uRy,u < 1, 

Condition (8) holds if and only if 



max min f(uB(w,UN),UN) 



w^B-^ -ulN^B-^\uJ, 



B-^w - B-^Nun 



< 1 



(9) 



When we consider the illustrative example in Section 1, we have B — [1], = [1] then 
problem (9) becomes 



max min f {ub{w , U2) , U2) = 

-1<W<1 U2&M. 



= - U2\U2] 



1 




W - U2 


1 




U2 



~ {w — U2Y -\-u\ <1 



(10) 



Now consider, function f{uB{w, mat), m^v)- It is a differentiable convex function in un- Then, for 

any w &W we can analytically determine the best response u*j^{d) — arg min^j^gigm-n f{uB{w, un),un), 
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by imposing 



Ru 



B-^w - B-^Nun 
un 



= 0, 



where / is the (m — n) x (m — n) identical matrix. We obtain 

-1 



u*n{w) = 



Ru 



-B~^N 
I 



-N^B-^^II 



Ru 



B-' 




w ~ —Mw, 



M 



where is the {m — n) x n null matrix. In the example under consideration, we have 

U*2{w) = 





1 




-1 




1 




1 




1 




1 




1 








w 



w 



For any w eW the minimal value of f{uj3{w, un),un) is 
where 



+ B-^NM 
-M 

V 

H 



H RuH 



(11) 



is a positive definite n x n matrix, as M is full rank. So far, we have shown that we can find 
the optimal value of problem (9) by solving problem 



max = w ^w, 

and checking that the optimal value is less than one. 



(12) 



We are ready to observe that problem (12) is easy as it reduces to determining the eigenvectors 
of an n X n matrix. 

Theorem 1 System (1) is e-stabilizable if and only if w*^^w* < 1, for all w* eigenvec- 
tors associated to the maximum eigenvalue of matrix R^^^ whose weighted quadratic norm 
w*^RwW* is equal to 1. 

Proof. As w^^w is convex, its optimal value w* lays on the frontier dW of the set W, i.e., for 
w*^RwW* — 1. Imposing the Karush Kuhn Tucker first order optimahty condition, we obtain 
2($ — XRyj)w* — 0. Then the optimal values of w* are some of the matrix R~^^ eigenvectors 
whose weighted quadratic norm w*'^RyjW* is equal to 1. In particular, w* are the eigenvectors 
associated to the maximal eigenvalues of R^^^- 
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□ 

In the example under consideration $ = | and w* — ±1 then w*^^w* = | < 1, hence the 
associated system is £-stabihzable. 

In the following we discuss for which initial state the system is certainly £-stabilizable through 
a (pure) linear state feedback control; hence we show that if we saturated the previous linear 
policy the system is £-stabilizable for any initial state. 



4 Ellipsoidal constraints 



Let us start by considering only the constraints (2) on w and neglect the ellipsoidal con- 
straints (3) on u. Among the saturated linear state feedback control (5) we prove that wc can 
solve Problem 1 using controls of type u = sat{—kHx}, with /c G M and if G M" as defined 
in (11). Note that matrix if is a right inverse of B, that is BH = I. We motivate the choice 
oiu — —sat{kHx} with H as defined in (11) as such a control describes the best response of 
u under the worst w as proved in the previous section. Also, note that the scalar A; G R must 
be lower than a certain value, which means that we cannot use a bang-bang control. This is 
motivated by the following reason. If we use a control u — sat{—kHx}, then the necessary 
and sufficient condition (8) becomes 

BUiin D W (13) 

where 

Klin = {m G : u = -kHx, k'^x^H^ R^Hx < 1}. 

Following the derivation of (12) in the previous Section, we have that (13) holds if and only 
if 

k'^w*^<^w* < 1. 

For k — 1 the above condition holds true as it reduces to (12). Obviously, the value k — 
\J w*'^^w* upper bound for k, namely, we must choose k such that A; < ^ if we wish the 
necessary and sufficient condition (13) be satisfied. 

With the above considerations in mind, we can conclude that the dimensions of the target 11 
where it is possible to drive the state are lower bounded. 

Denote by XmaxiZ) the maximum eigenvalue of a given matrix Z. In the following theorem 
we prove that V{x) < within a given set {invariant set). This result will allow exploiting 
V{x) as a Lyapunov function to prove the convergence to the target set 11. 

Theorem 2 Consider system (1) subject to the only ellipsoidal constraints (2) on w, and 
controlled via linear state feedback u — —kHx, with H such that BH — I. Then condition 
V <Q holds if and only if 

k'^{x^Pxf - x^PR-^Px > 0. (14) 
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Proof. For H such that BH — I, condition 1/ < is equivalent to 

2kx^Px + 2w^Px > 0. (15) 

We aim at proving that V < holds for any x external to an appropriate smooth closed 
surface. To do this, we look for an x e R" inducing a solution strictly greater than zero for 
the following problem 

min ((x, w) = 2kx^Px + 2w'^Px. (16) 

As C(x, w) is linear in w, the optimal w* must lay on the boundary of set W. The Karush Kuhn 
Tucker conditions impose that Px = —XRwW* for some A > 0, that is w* = —jR~^Px. Note 
that being P full rank, it necessarily holts thatA 7^ for all x 0. Then, ({x, w*) = 2kx'^Px — 
jx'^PR'^Px > 0. As w* lays on the boundary of W, we have w*'^RyjW* = ^'^pR^ ^ 1 

from which A = ^Jx'^PR~^Px. Hence, C,{x, w*) > 0, and therefore also (15) holds, if and only 
if (14) holds. 

□ 

We now exploit V{x) = x^ Px as a Lyapunov function to prove the convergence to the target 
set n. We determine under which conditions on P and k we have that y < or, equivalently, 
inequahty (14) hold for any x ^ 11. 

When P = uRw, (14) becomes k'^x^Px > u. Then, in this case, we can use V{x) to prove 
the convergence of the system to 11 for /c^ > u. 

In the following, we consider the general case when P ^ i/Ryj. 

Lemma 1 Consider system (1) subject to the only ellipsoidal constraints (2) on w, and 
controlled via linear state feedback u = —kHx, with H such that BH = I . Then, k'^{x'^ PxY — 
x^PR~^Px > holds for any x if and only ifk^ — x^PR~^Px > holds for any x e 911. 

Proof. (Necessity) . Assume that there exists x G dU such that k"^ — x^PR^^Px < 0. Then, 
there also exists a ball Ball{x,r) centered in x with a sufficiently small radius r > such 
that for all x G Ball{x, r) we have k^ — x^ PR^^Px < 0. This implies that there exist a; ^ 11 
for which condition (14) does not hold. 

[Sufficiency) . Assume that k'^ — x^PR^^Px > holds for any x G dU. By contradiction, 
consider x ^ IT, i.e., x^Px = p > 1, such that k'^{x^Px)'^ — x^ PR^^Px < 0, that is 
k'^p'^ - x^ PR:^^Px < 0. Then, there exists x = ^ e dU such that /c^p^ - px^PR:^^Px < 0, 

that is k'^p — SFPR~^Px < 0. This latter result is contradictory as we cannot have k'^p < 
x^PR-^Px < A;^ for p > I. 

□ 

Lemma 2 Consider system (1) subject to the only ellipsoidal constraints (2) on w, and 
controlled via linear state feedback u — —kHx, with H such that BH — I. We can use V{x) 
to prove the convergence of the system to U for k^ > X^axiRw^P)- 

Proof. Condition k^ — x^ PR~^Px > holds for any x G dU if and only if min^g^nj^^ — 
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x^PR'^Px} > 0. Imposing the Karush Kuhn Tucker first order optimality condition, we 
obtain 2(PR~^P — XP)x* = 0. Then the optimal vahics of x* are some of the matrix R^^P 
eigenvectors whose weighted quadratic norm x*^Px* is equal to 1. In particular, x* are 
the eigenvectors associated to the maximal eigenvalues of R~^P- For vectors x*, condition 
A;2 - x*^PR-^Px* > becomes k'^ - Xmax{K^P)x*^Px* > 0, that is k'^ - XmaxiK^P) > 0- 

□ 



Observe that the system converges to the target set Ilij = {x : k'^x^ Ru,x < 1} as any feasible 
target set H = {x : x^Px < 1}, with k"^ > Xmax{Rw^P) includes Ur. Indeed, 11 ^ 11^ if 
x'^Px — k'^x'^RyjX — x'^{P — k'^Ryj)x < or equivalently if P — k'^Ryj :< 0. In turn, the latter 
condition is equivalent to R~^P — k"^! :< that certainly holds as k'^ > X^axiRw^P) 

In the next theorem we introduce the constraints on controls (3). To this end, we need to 
define the family of ellipsoid So(^) = e M" : x^Px < x{0)^Px{0) := ^} parametrized in 

Theorem 3 Given system (1) in the ellipsoidal case, we can drive the state x{t) from any 
initial value x{Q) G So(^) to the target set 11 via linear state feedback u — —kHx if the 
following conditions hold 

> X^ax{R'JP) (17) 

k-'iX^axiP-'^) < 1. (18) 

Proof. By Lemma 2, under condition (17) it holds V{t) < for all x{t) ^ U and then V{x) 
can be considered as a Lyapunov function for the convergence of the state to the set 11 when 
the linear control u = —kHx is implemented. Condition V{t) < also implies that ^o{C,) is 
invariant with respect to the same linear feedback as ^ > 1 which means So(0 3 11- Then 

maxu^(t)Ruu(t) < max k^x'^ H'^ RuHx = max k'^x'^^x ^ k^^Xmax(P~^^)- 

t>o xeT.o{0 xeT.o(0 

Therefore the constraint u — —kHx{t) e U for alH > is enforced if (18) holds true. 

□ 



The following theorem provides a solution to Problem 1. Let us denote by X the set of states 
X where we can define a linear control u{x) — —kHx, i.e., X — {x : —kHx e U}. Consider 
the saturated hnear state feedback control of type 



u{x) 



-kHx 



Hx 



if a; G X 
iix^X 



(19) 



Theorem 4 Consider a system (1) in the ellipsoidal case. For any positive definite matrix 
P e R"^" satisfying condition (17), the saturated linear state feedback control (19) drives the 
state x{t) within the target set U for any initial state x{0). 

Proof. By construction, u{x) is a continuous function with U as codmain. When we use such 
a control, we know that V{x) < also holds for any x ^ 11, if 11 C X and k'^ > Xmax{R~^P) 
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(see Lemma 2). 



First observe that, for all x G dX, we have x^Px > k'^x^H^RuHx = 1, where the lat- 
ter inequality holds as 11 C X. Then, for any x ^ X, that is for k^x^H^RuHx > 1, we 
x^htr^Hx > ^ ^max{R~'^P) siuce both x'^Px and x^H'^RuHx are positive definite 
quadratic forms. 

In Lemma 2, we have proved that V{x) < for x G X \ 11. Now, we consider x ^ X. We 
have V{x) < if and only if —x^PBu{x) + x^Pw > 0, for all w G W, that is 



Px 1 
, ,„ + x'^Pw } > (20) 

must hold. Applying the Karush-Kuhn-Tuker conditions, we transform (20) in x'^Px _ 

Y x^ H"^ RuHx 

^Jx'^P'^R^'^Px > 0. In turn, the latter inequality holds if ^thtr^hx ~ ^maxiR~^P) > 0, as 
x'^P^R-^Px < Xmax{R''^P)x^Px. Wc then conclude that V{x) < since ^thIr^hx > > 

□ 

Observe that the saturated linear state feedback control (19) is not decentralized in the sense 
that the generic ith control Ui in general depends on the demand at different nodes and on 
the other controls Uj, j ^ i. This is due to either the structure of matrix H or the ellipsoidal 
constraints (3). 

Remark 1 Consider the two equivalent matrix inequalities on P and Q = P~^ , 

{2k - 1)P - PR-^P > 0, {2k - 1)Q - R-^ < 0. (21) 

Trivially, any P satisfying condition (17) also satisfies the two above matrix inequalities. 

Matrix inequalities of the above form will be used in the following sections. 

Example 2 Consider the graph depicted in Fig. 1, with one node and two arcs and incidence 
matrix B = [1 1]. Controls are subject to ellipsoidal constraints (6). Then we have, R^ = 1, 
Ru — I and $ = i. We can stabilize the system within U — {x & M. : x^ < 1} for any initial 
state x{0) < V2 via a pure linear state feedback u — —kHx. To see this take Q — I, and 
observe that the matrix inequality on Q (21) is satisfied for any k > 1. Furthermore, if we 
assume k — 1, then from (18) we must have k"^ = 1 < = 



5 Polytopic constraints 



Controls u are subject to the polytopic constraints (4). Again, we study under which condi- 
tions we can solve Problem 1 using controls of type u — —sat{kHx}, with /c G M and i7 G M** 
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such that BH — I. In this case, we interpret the sat{.} operator as componentwise. More 
specifically, we choose the control 



Ui 



sat[^-,u+]{-kHi,x}, 



(22) 



with H such that BH — I, Hi, denoting the ith row of H and where, for any given scalar a 
and b 

' 6, a c>b, 

C, if a < C < 
a, if C < 

Henceforth we omit the indices of the sat function. 



saha,b]{C} = ' 



Under the control u — sat{—kHx}, the closed loop dynamics becomes 

X = Bsat{-kHx} - w. (23) 
Our idea is to rewrite the above dynamics in the following polytopic form 

X = A{t)x{t) - w{t), w{t)'^R^w{t) < 1, (24) 

where the time varying matrices A{t) are expressed as convex combinations of 2"^ matrices 
Aj, j = 1, . . . ,2"^. More precisely the expressions for A{t) are 



(25) 



The procedure to compute matrices Aj's is borrowed from [9] and recalled below. Let us 
rewrite the control policy as 

Ui — sat{—kHi,x} = 9i{x){—kHi,x), 

where 9i{x) are the "degree of saturation" of the control components defined as follows 

if —kHi,x < uj 



Oi{x) 



1 ii < —kHi,x < uf 
— r# — if —kHi,x > u'^ 



(26) 



Let 6_ = [6_i, . . . ,£„J be a vector whose components arc such that < £j < 1 and represent 
lower bounds of 6i{x{t)), for t > 0. Lower bounds depend on a;(0) and can be computed 
as 6_i = minj.gx;o(0 ^i(^) where we remind the definition of ^o{C,) = {x E M" : x^Px < 
x{0)^Px{0) := Also define the vector ip^ = [ip^, . . . ,ip^] with ipf = j- and the associated 

—i 

portion of the state space 



S{iP^) = {x e M" : -V^^ < -kHx < 
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According to the above definition of the 9iS we derive that 2 5]o(0- Note that we can 

affirm that 9^ are lower bounds because the state trajectory never exits S{ip^) as we will show 
in the proof of Theorem 5. 

Consider now the 2"* vectors 7j G {1, fii} x . . . x {1,9^}, with j = 1, . . . , 2™. In other words, 
7j is an m component vector with ith component taking value 1 or ^j. Then, each matrix 
Aj can be expressed as Aj — —Bkdiag{jj)H. Roughly speaking each vector 7, stores the 
minimum and or maximum degree of saturation of all control components. Also, note that 
matrices AjS induce a partition of S{ip^) into regions Xj, with j = 1 . . . , 2"\ Each region is 
defined as the set of state values such that the control components are saturated with degree 
of saturation equal to 7jj, namely 

Xj ^ {x eW^ : 9i{x) = 7ji, i = 1, . . . , m}. 

We remind here that 7^1 is the ith component of 7^ . 

To complete the derivation of the polytopic form (24) it is left to be noted that given any 
x{t) G S{iIj^) we can compute the associated degree of saturation from (26) and derive the 
weights crj{t) of the convex combination (25). All the results in the rest of this section try to 
give an answer to Problem 1 with respect to the polytopic system (24). For each Aj, let us 
define a matrix 

Mj = QAj + AjQ + aQ + -R'^ 

for a given positive and arbitrarily chosen scalar a and let (A^-, wj) with re {1, . . . , n} be the 
negative eigenvalues and corresponding eigenvectors of Mj. 

Theorem 5 Consider system (1) in the polytopic case. The saturated linear state feedback 
control (22) drives the state x{t) within the target set H if 

Xj C Span{v'j}, for all j = 1, . . . , 2". (27) 



Proof. First of all, note that if (27) holds true then ^o{C,) is invariant. Consequently, as 
So(0 ^ S{ip^) and by definition x{0) G So(0' ^^^^ have that the state trajectory x{t) 
will never exit S{'ijj^). Now, we must show that V{x) < for all x and w such that x ^ U, 
u & U and w G W. In formulas, we must have 

V{x) = x^Px + x'^Px = [A{t)x - wfPx + x'^P[A{t)x - w] ^ ^^^^ 
= x'^A{t)'^Px + x'^PA{t)x - w^Px - x'^Pw < 

for all X and w satisfying 

1 - x'^Px < (29) 
w'^RyjW - 1 < 0. (30) 

Using the »S-procedure, we can say that condition (28) is implied by conditions (29)- (30) if 
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there exist a, /3 > 0, such that for all x and w 



X 


T |- 


w 





A(tfP + PA(tf + aP -P 
-P -PRu 





X 




w 



a + /3<0. 



(31) 



Trivially it must hold (3 < a. Assume without loss of generality (3 = a. Remind that a and 
(3 can be chosen arbitrarily. After pre and post-multiplying hy Q = P~^, the above condition 
becomes ^ 

'QA{tf + A{tfQ + aQ -I 



X 


T [- 


w 





-I 



—aRu 





X 




w 



< 0. 



(32) 



Now, as the state never leaves the region S{ip^), i.e., x{t) G S{ip^), we can always express 
A{t) as convex combination of the AjS as in (25). 



By convexity, the above condition is true if it holds, for all j = 1, . . . , 2", 





T r 


^(i) 









QAj + AjQ + aQ -I 
—I —aR,, 



< 0. 



(33) 



Using the Shur complement the condition (33) is implied by (27). 

□ 

Stronger conditions arc established in the following theorem which also highlights the depen- 
dence of Mj on the scalar a. 

Theorem 6 Consider system (1) in the polytopic case. The saturated linear state feedback 
control (22) drives the state x{t) within the target set U if there exists a scalar a > such 
that 

Mj<0, forallj = l,...,2'\ (34) 
Proof. Trivially, if we observe that (34) implies (27). 

□ 

Both (27) and (34) are sufficient, but not necessary, conditions. When they hold, we are sure 
that the system state converge to a state strictly included in the target set 11. We discuss 
more on this topic in the next section. 



5.1 Approximation error 



Wc wish to estimate the difference in terms of volumes between the target set 11 and the target 
set obtained from conditions (34) and we will call such a difference as approximation error. 
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On this purpose, denote by Qj the matrix of the smallest (in volume) ellipsoid satisfying 
Mj < 0, which is given by 

Qj = arginf min{ciei(g), Mj = QAj + AjQ + aQ + -R'^ < 0}. (35) 

To do this, let matrix A be the matrix Aj with j = 1, . . . , 2"^ obtained when no controls are 
saturated and note that the dynamics associated to this single matrix is the same as if we 
assumed the controls unbounded. To be more precise, A — —BkH as all components of 7^ 
are equal to one. Remind that 7j stores the degree of saturation of each control component. 
Also let us define Q the solution of (35) for Aj = A. Wc do this, as the target set 11 within 
which we can stabilize the state, must inscribe the ellipsoid defined by Q, i.e., 

n D {x e R" : x'^Q-^x < 1}. 

Similarly, let matrix A be the matrix Aj with j = 1, . . . , 2™ obtained when all controls are sat- 
urated at their lowest degree of saturation. To be more precise, A = —Bkdiag{[6i, . . . , 9jn])H 
as all components of 7^ are equal to for i = 1, . . . ,m. If we also define Q the solution 
of (35) for Aj — A, the target set 11 must be inscribed in the ellipsoid defined by Q, namely, 

n C {x e M'' : x'^Q~\ < 1}. 



The approximation error can be measured by the ratio 

det(Q~^) - det{Q~ 



e — 



det{Q-^) 



Example 3 Consider the graph depicted in Fig. 1, with one node and two arcs, incidence 



matrix B = [1 1], and target set 11 = {x G 



: X 



constraints (1). Take H = [ 
scalar) 



1 llT 



< 1}. Controls are subject to polytopic 



and k — 1. Then according to (26) we have (here x is a 



Oi{x) 



^/^/2>2 
1 if-3<x/2<2 
-^,ifx/2<-3 



02{x) 



1 if -l<x/2< 2 
if x/2<-l 



If we consider initial states x{0) satisfying — 10 < a;(0) < 10, possible lower bounds for the 
0's are 9i = I and 02 = ^. Note that S{ijj^) = {x e : -10 < x < 10}. Vectors 7^5 and 
matrices A 's turn out to be 



71 = [1 1^72 = [0.4 ir73 = [l 0.2]^ 74 = [0.4 0.2f 
Ai = -2 A2^ -1.4 A3 = -1.2 A4 = -0.6 

Dynamics (24) is then 

X = [-ai{t)2 - (72(^)1.4 - (73(^)1.2 - a4{t)0.6]x + w, 



(36) 



(37) 
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with Yl^=i(Tj{t) — 1. Furthermore, we have 

Mi = [-4 + Q;]g + ^ M2 ^ [-2.8 + a]Q + ^ 
M3 = [-2.4 + a\Q + lM,^ [-1.2 + ajQ + ^ ' 

To apply Theorem 5 and 6, note that A = A^ and that M4 < implies consequently Mj < 
for all j. The solution of (35), forj = 4 Z5 Q4 = Q = and a = 0.6, then the approximation 
error is e = = 1.78. ~ 

6 Conclusions and future works 

We have addressed the problem of e-stabihzing the inventory of a continuous time hn- 
ear multi-inventory system with unknown demands bounded within ellipsoids and controls 
bounded within ellipsoids or polytopes. Motivations are due to the cost reduction associated 
with a bounded inventory. As main results we have provided certain LMIs conditions under 
which £-stabilizability is possible through a saturated linear state feedback control. We have 
also exploited some recent techniques for the modeling and analysis of polytopic systems with 
saturations. 

This work is a continuation of [2] and is in line with some recent applications of LMI tech- 
niques to inventory/manufacturing systems [6]. In a future work, we will study the validity 
in probability of the LMI conditions derived in this paper. This is in accordance with some 
recent literature on chance LMI constraints developed in the area of robust optimization [4,8]. 
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